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L EXIN  GTON 


M V S S A C H U SETTS 


ABSTRACT 


As  part  of  the  work  on  the  development  of  a submarine-towed  ELF  loop 
antenna,  the  mechanics  of  cable  vibration  have  been  studied  both  analytically 
and  experimentally.  In  particular,  measurements  of  towed  cable  vibrations 
have  been  carried  out  from  surface  vessels  and  from  a submarine.  This  note 
presents  the  results  of  these  measurements,  which,  after  some  preliminary 
accelerometer  measurements,  were  obtained  with  a buoyant  cable  equipped  with 
strain  gages. 

The  results  are  that  both  cable  longitudinal  strain  and  cable  curvature, 
over  the  frequency  range  of  20  to  200  Hz,  are  excited  by  the  fluctuating 
forces  on  the  surface  of  the  cable  arising  from  the  cable’s  own  turbulent 
boundary  layer.  Neither  vibration  of  the  tow  point  nor  the  submarine  wake 
are  significant  sources  of  mechanical  energy.  Generalized  force  spectra  are 
derived  from  the  measured  vibration  spectra.  These  then  allow  the  levels  of 
both  types  of  vibration  to  be  calculated  for  an  arbitrary  cable  towed  at  an 
arbitrary  speed. 
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Strain-Gage  Vibration  Measurements  on  a 
Submarine-Towed  Antenna  Cable 


I.  Introduction 


This  report  describes  the  results  of  strain-gage  measurements,  made  in 
1972,  of  the  vibration  of  a submarine-towed  antenna  cable.  These  measurements, 
taken  together  with  some  preliminary  measurements  made  on  a cable  fitted  with 
an  accelerometer  array,  have  played  a decisive  role  in  the  analysis  of  towed 
antenna  noise  and  in  the  development  of  techniques  to  overcome  it  [1,2, 3, 4, 5]. 

The  significance  of  knowing  how  and  why  the  cable  vibrates  can  be  illus- 
trated by  the  following  questions: 

i)  How  long  should  the  tail  be?  (The  tail  is  the  electrically 
inactive  length  of  cable,  aft  of  the  antenna  section,  and  used  to 
isolate  the  antenna  from  any  end  instabilities.) 

ii)  Will  it  help  to  place  mechanical  isolators  in  the  cable  to 
prevent  vibrational  energy  propagating  along  the  cable  from  tail  or 
tow  point  into  the  antenna  region? 

iii)  Would  the  use  of  a longer  antenna  cable  help?  The  longer 
cable  would  place  the  antenna  region  of  the  cable  further  away  from 
the  large  wake  eddies  of  the  submarine. 

iv)  Is  the  cable  vibration  the  only  significant  source  of  the 
noise  induced  by  towing  the  antenna,  or  is  hydromagnet ic  noise  a 
factor?  Hydromagnetic  noise  is  the  noise  voltage  induced  in  the 
antenna  by  the  electrical  currents  in  the  surrounding  water.  These 
currents,  in  turn,  are  produced  by  the  turbulent  motion  of  the 
electrically  conductive  sea  water  in  the  geomagnetic  field. 

The  next  section  describes  briefly  the  results  of  accelerometer  measure- 
ments, which  provided  invaluable  experimental  support  in  the  evolution  of  the 
theory  of  vibrational  wave  propagation  along  the  towed  cable.  However,  the 
accelerometer  installation  unavoidably  so  modified  the  cable  diameter  and 
mechanical  properties  that  it  could  not  pass  out  through  the  deployment  mech- 
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anism  of  a submarine  and  would  not,  in  any  case,  give  an  accurate  measure  of 
the  vibration  level  of  an  unmodified  cable.  The  measurements  were  therefore 
carried  out  from  a surface  boat  using  a depressor  (a  device  like  an  upside 
down  kite)  to  lower  the  tow  point  to  the  appropriate  depth. 

The  basic  quantitative  vibration  data  for  an  unmodified  cable  towed  behind 
a submarine  were  obtained  with  strain  gages.  These  measurements  are  described 
in  Section  III  and  the  results  are  digested  and  summarized  in  Section  IV. 

II . Accelerometer  Measurements 

The  accelerometer  experiment  was  carried  out  as  a joint  venture  by  the 
Navy  Underwater  Systems  Center,  New  London  Laboratory  and  Lincoln  Laboratory. 
The  latter  designed  the  experiment  and  directed  the  data  processing.  The 
former  was  responsible  for  installing  the  accelerometers,  for  providing  the 
instrumentation  and  for  carrying  out  the  data  processing. 

A total  of  seven  accelerometers  were  installed,  an  array  of  five  to  mea- 
sure transverse  accelerations  in  a common  plane  at  five  different  points 

within  a total  distance  of  3.4m,  one  to  measure  the  transverse  acceleration 
in  the  perpendicular  direction  and  one  to  measure  in  the  longitudinal  direc- 
tion (see  Fig.  1).  The  towing  configuration  is  shown  in  Fig.  2. 

The  accelerometer  signals  were  recorded  on  magnetic  tape,  digitized  and 
then  processed  in  pairs  to  obtain  the  individual  auto  spectra,  the  cross 
spectrum,  the  coherence  function,  the  transfer  function  and  the  impulse 
response . 

In  the  first  experimental  run,  the  most  striking  feature  in  the  result 
of  processing  between  pairs  in  the  array  of  five  was  the  disclosure  that  the 
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Fig.  1.  The  placing  and  orientation  of  the  seven  accelerometers  in  the  cable. 
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Fig.  2.  The  towing  configuration  of  the  boat,  depressor  and  cable. 
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transverse  waves  were  propagating  in  one  direction  only — towards  the  tow  point. 
Figure  3 shows,  for  example,  the  impulse  response  between  the  two  accelerom- 
eters 3. Am  apart.  If  the  cable  vibration  had  had  transverse  waves  of  compar- 
able amplitude  traveling  in  both  directions,  there  would  be  a structure  on 
the  right  side  of  the  time  origin  mirroring  that  developed  very  clearly  on 
the  left. 

The  shape  of  the  impulse  response  is  that  expected  for  a dispersive 
transmission  line.  The  wave  velocity  and  attenuation  are  greater  at  higher 
frequencies,  due  to  the  effect  of  cable  stiffness,  so  that  the  smaller  ampli- 
tude higher  frequency  waves  arrive  first. 

The  phase  response  is  shown  in  Fig.  4,  and  this  again  gives  strong  evi- 
dence of  the  unidirectional  propagation,  for  the  phase  is  a monotonically 
increasing  function  of  frequency.  It  can  be  used  as  a check  on  the  stiff  taut 
string  model  of  the  mechanical  properties  of  the  cable,  which  expresses  the 
relationship  between  the  wave  number  of  the  propagating  wave  and  the  ratios 
stiffness  over  mass  and  tension  over  mass  [2].  In  Fig.  5,  the  plotted  points 
were  obtained  by  noting  the  frequencies  at  which  successive  half-cycle  incre- 
ments in  phase  occur  in  Fig.  4.  The  continuous  line  is  the  phase  difference 
predicted  by  the  stiff  taut  string  model  using  the  parameter  values  shown. 

Thus  the  general  form  of  the  propagation  seems  to  be  well  described  by  the 
stiff  taut  string  model. 

It  should  be  noted  that  the  accelerometer  installation  caused  a very 
marked  stiffening  of  the  cable  and  that  the  tail  end  of  the  cable  tangled  into 
a large  knot,  which  was  not  noticed  until  the  cable  was  being  retrieved,  and 
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Fig.  3.  The  impulse  response  between  the  two  accelerometers  3.4  m apart  in 
the  array  of  five. 
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Fig.  4.  The  phase  response  between  the  two  accelerometers  3.4  m apart  in 
the  array  of  five. 
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Fig.  5.  Phase  difference,  for  transverse  wave  propagation  over  a distance 
of  3.4  m,  versus  frequency. 
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which  would  make  the  tension  larger  than  expected.  It  is  not  surprising, 

therefore,  that  the  parameter  values  used  to  make  the  curve  fit  the  measured 

data  are  larger  than  those  obtained  from  direct  measurement  of  the  physical 

A “ 2 39  _2 

properties  of  the  cable,  which  are  3.4  m sec  ^ for  El/m  and  0.62  x 10  m sec 
for  T/m. 

Other  results  of  this  first  run  were  that  the  coherence  between  mutually 
orthogonal  pairs  of  accelerometers  was  in  all  cases  too  small  to  measure, 
indicating  that  the  vibrations  in  the  three  principal  directions  are  statis- 
tically independent.  In  addition,  the  output  levels  of  all  transverse  sen- 
sors remained  constant  in  time,  which  indicates  that  the  transverse  vibration 
was  statistically  isotropic.  For  if  it  were  not,  then  the  sensors  would  have 
shown  a signal  level  variation  due  to  cable  rotation. 

For  this  first  test,  the  300  ft.  tail  section  was  mechanically  fastened 
to  the  rest  of  the  cable  with  a braided  metal  sleeve  clamped  to  each  cable 
piece  with  two  standard  hose  clamps.  The  clamps  were  located  about  five  feet 
aft  of  the  last  sensor.  It  was  suspected  that  the  protruding  tightening 
screw  on  each  clamp,  which  generated  hydrodynamic  disturbances  in  the  water 
and  was  in  turn  buffeted  by  them,  was  the  source  of  the  cable  vibrations. 

Their  direction  of  energy  flow  was  certainly  consistent  with  this  view. 

For  subsequent  tests,  therefore,  the  tail  was  reconnected  using  a smooth- 
ly faired  connector.  The  results  were  i)  the  level  of  the  vibration  dropped, 
ii)  the  phase  and  cross-spectrum  taken  between  pairs  of  accelerometers  showed 
the  features  characteristic  of  essentially  equal  amplitude  uncorrelated 
energy  flow  in  both  directions,  iii)  shortening  the  tail  by  100  ft.  caused  a 
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slight  reduction  in  the  level  of  the  vibration  but  left  the  symmetry  of  the 
energy  flow  unchanged.  Thus,  in  view  of  the  large  attenuation  of  transverse 
waves  [2],  the  mechanical  energy  must  arise  locally,  for  otherwise  the  sym- 
metry of  the  energy  flow  would  have  been  profoundly  changed  by  what  was,  in 
effect,  a change  in  the  location  of  the  array  of  100  ft.  In  addition,  the 
view  that  the  hose  clamps  were  the  source  of  the  vibration  measured  in  the 
first  test  was  shown  to  be  true. 

The  phase  plots  between  the  same  two  accelerometers  for  the  first  run  and 
for  the  second  run  at  both  tail  lengths  are  shown  in  Figs.  6,  7 and  8. 

Clearly  marked  is  the  difference  between  the  smooth  phase  progression  of  the 
first  and  the  progression  by  jumps  of  the  other  two.  The  corresponding 
cross-spectra  are  shown  in  Figs.  9,  10  and  11.  Again  the  difference  between 
the  first,  which  show  the  characteristic  periodicity,  is  clearly  marked.  The 
minima  in  the  cross-spectra  coincide  with  the  phase  jumps  in  the  correspond- 
ing phase  plots,  as  expected. 

No  longitudinal  acceleration  data  are  reported  here  because  the  attenua- 
tion of  longitudinal  waves  in  the  cable  is  small  [2].  This  allows  disturb- 
ances generated  at  the  depressor  to  travel  along  to  the  accelerometer  location 
and  dominate  other  sources  of  longitudinal  strain  energy  there.  Also,  since 
only  one  longitudinally  oriented  accelerometer  was  used,  no  longitudinal  wave 
structure  analysis,  like  that  carried  out  for  transverse  waves  between  pairs 
of  transversely  oriented  accelerometers,  was  possible. 

The  outcome  of  the  accelerometer  measurements  was,  then,  that  the  two 
orthogonal  directions  of  the  transverse  vibration  and  the  longitudinal  vibra- 


10 


Fig.  6.  Phase  response  at  2.2  m spacing  and  8 knots  with  a 300  ft.  tail  and 
hose  clamps  present. 
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Fig.  7.  Phase  response  at  2.2  m spacing  and  8 knots  with  a 300  ft.  tail  and 
no  hose  clamps  present. 


12 


Fig.  8.  Phase  response  at  2.2  m spacing  and  8 knots  with  a 200  ft.  tail  and 
no  hose  clamps  present. 
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Fig.  9.  Cross-spectrum  at  2.2  m spacing  and  8 knots  with  a 300  ft.  tail  and 
hose  clamps  present. 
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CROSS-SPECTRAL  DENSITY  ( dB,  arbitrary  origin) 


Fig.  10.  Cross-sprectrum  at  2.2  m spacing  and  8 knots  with  a 300  ft.  ta 
and  no  hose  clamps  present. 
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Fig.  11.  Cross-spectrum  at  2.2  m spacing  and  8 knots  with  a 200  ft.  tail 
and  no  hose  clamps  present. 
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tion  are  all  statistically  independent  from  one  another,  that  the  propagation 
of  transverse  vibrational  waves  along  the  cable  is  well  described  by  the  stiff 
taut  string  mathematical  model  and  that  the  transverse  vibration  originates 
neither  at  the  tail  end  nor  at  the  tow  end  but  is  locally  generated.  Whether 
the  boat  wake  or  the  turbulent  boundary  layer  of  the  cable  itself  provides 
the  disturbing  transverse  force  cannot  be  decided  from  these  measurements. 

Ill . Strain-Gage  Measurements 

The  strain  gage  installation  was  intended  to  provide  a means  for  measur- 
ing the  longitudinal  and  transverse  vibration  of  the  cable  separately  and  at 
the  same  time  leave  the  mechanical  properties  and  dimensions  of  the  cable 
unchanged.  The  installation  and  instrumentation  have  already  been  described 
by  Nackoney  [6].  Briefly,  there  are  six  gages  installed  as  three  pairs  at 
three  locations  in  a 2000  ft.  long,  14  conductor  buoyant  cable  (Fig.  12a). 

The  two  gages  in  each  pair  are  emplaced  diametrically  opposite  one  another  in 
the  buoyant  cable  jacket  and  have  their  sensitive  axes  parallel  to  the  cable 
axis  (Fig.  12b).  The  sum  of  the  strains  registered  by  each  member  of  the 
gage  pair  is  a measure  of  the  longitudinal  strain  of  the  cable  at  the  gage- 
pair  location.  The  difference  between  the  strains  is  a measure  of  the  curv- 
ature of  the  cable  axis  at  the  same  location. 

If  the  gages  are  perfectly  matched  and  the  installation  is  both  geometric- 
ally and  mechanically  perfectly  symmetrical,  the  strain  sum  should  be  independ- 
ent of  the  cable  curvature  and  the  strain  difference  should  be  independent  of 
the  cable  longitudinal  strain.  But  the  gages  are  not  perfectly  matched,  the 
cable  is  not  perfectly  symmetrical  and  the  installation  procedure  is  imper- 
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Fig.  12.  Strain  gage  experiment  showing  gage-pair  location  and  the  gage 
placement  in  the  cable  cross  section. 
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feet.  So  the  question  arises  of  how  the  gain  of  the  amplifier  serving  each 
member  of  a gage  pair  is  to  be  determined  to  reject  the  unwanted  signal.  This 
question  is  particularly  critical  if  the  curvature  signal,  for  example,  is 
very  much  larger  than  the  longitudinal  strain  signal.  For  then,  a slight 
error  in  gain  setting  in  forming  the  strain  sum  could  lead  to  enough  of  the 
curvature  signal  left  uncancelled  to  obliterate  the  required  longitudinal 
strain  signal.  It  is  also  clear,  that,  due  to  differences  in  gain  and  phase 
of  the  two  amplifiers,  the  proper  gain  setting  could  well  be  a complex  number 
which  is  a function  of  frequency. 

Fortunately,  the  longitudinal  and  transverse  vibrations  of  the  cable  can 
be  assumed  to  be  statistically  independent.  Such  an  assumption  was  shown 
experimentally  to  be  valid  for  the  accelerometer  array  cable.  It  is  shown 
in  the  Appendix  that  this  simple  assumption  makes  it  possible  to  define  a 
method  of  data  processing  which  establishes  the  proper  gain  setting,  as  a 
function  of  frequency,  in  phase  and  amplitude,  by  means  of  certain  averages 
performed  on  the  two  strain  signals.  This  gain  adjustment,  and  its  applica- 
tion to  calculate  the  separate  longitudinal  strain  and  curvature  spectra,  was 
carried  out  in  a digital  computer.  The  input  to  the  computer  was  the  two 
separate  digitized  and  synchronized  strain  signals.  Its  output  was  the  sep- 
arate required  spectra. 

The  first  test  of  the  strain  gage  cable  was  carried  out  on  the  surface. 

It  served  as  a dress  rehearsal  for  the  submarine  test  to  follow.  No  depres- 
sor was  used,  so  the  cable  simply  floated  along  the  surface  behind  the  towing 
boat.  The  test  was  successful  in  demonstrating  that  the  whole  system — cable, 
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instrumentation  and  recording — functioned  properly.  It  was  uneventful  except 
for  a valuable  corroboration  of  the  theory  for  longitudinal  strain  wave  prop- 
agation along  the  cable. 

Figure  13  shows  short  (c.  4.5  sec)  simultaneous  samples  of  the  sum  strain 
from  each  gage  pair.  These  were  obtained  by  using  summing  amplifiers  rather 
than  the  digital  process  described  in  the  Appendix  (which  can  produce  only 
spectral  averages).  Each  amplifier’s  gain  balance  setting  was  determined 
solely  by  the  criterion  that  the  gain  in  all  channels  be  the  same. 

The  strain  signals  show  very  strong  correlation  from  pair  to  pair  with  a 
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delay  equivalent  to  a speed  of  about  6 x 10  ft. /sec  apparent  between  the  tow- 
point  pair  and  the  tail  mounted  pairs.  The  tail  pair  signals  are  of  smaller 
magnitude  than  the  tow-point  ones.  Thus,  for  the  dominant  frequency  compon- 
ents visible  here  and  for  this  surface  test  of  the  cable,  one  concludes  that 
the  longitudinal  strain  energy  was  generated  near  the  tow  point  and  then  prop- 
agated towards  the  tail. 

Figure  14  shows  the  corresponding  simultaneous  samples  of  the  difference 
strain  from  each  gage  pair.  These  signals  are  proportional  to  the  local 
curvature  of  the  cable.  There  is  no  visually  apparent  correlation  from  pair 
to  pair  and  the  curvature  variation  increases  in  magnitude  towards  the  tail. 
This  suggests  that  the  curvature  (and  therefore  also  the  transverse  vibration) 
was  locally  generated,  for  this  surface  test. 

The  submarine  test  consisted  of  towing  the  strain  gage  cable  at  various 
speeds  in  straight  runs  for  three  different  cable  configurations,  executing 
a continuous  heading  change  (i.e.,  a circular  turn)  at  various  speeds  and 
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Fig.  13.  Simultaneous  longitudinal  strain  signals  from  the  three  gage  pairs 
at  a speed  of  8 knots  on  the  surface. 
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Fig.  14  Simultaneous  curvature  signals  for  the  three  gage  pairs  at  a speed 
of  8 knots  on  the  surface. 
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performing  a slow  speed-increase  run.  For  each  run,  the  strain  signals  from 
the  six  gages  were  recorded  separately  on  six  channels  of  a magnetic  tape 
recorder.  After  the  test,  the  signals  were  digitized  and  processed  as  des- 
cribed in  the  Appendix. 

The  three  different  cable  configurations  used  in  the  straight  runs  were, 
first,  with  1850  ft.  of  cable  being  towed  as  shown  in  Fig.  12;  second,  with 
only  1150  ft.  of  cable  being  towed  and,  third,  after  200  ft.  had  been  cut  off 
the  tail  end,  with  1650  ft.  of  cable  being  towed.  For  the  turns  and  for  the 
slow  speed  increase  runs,  the  cable  configuration  used  was  the  first  (1850  ft. 
total  tow  length,  300  ft.  from  tail  end  to  gages  5 and  6). 

The  story  told  by  the  longitudinal  strain  results  was  that  the  longitud- 
inal strain  spectrum  is  largely  the  same  at  all  gage  pair  stations  for  both 
the  long  scope  straight  runs  (1650  ft.  and  1850  ft.  of  cable  towed)  and  the 
short  scope  straight  run  (1150  ft.  of  cable  towed).  Figure  15  shows  the 
longitudinal  strain  spectra  at  12  knots  from  a total  of  eight  gage  pair  loca- 
tions superimposed  on  one  another.  The  locations  were  1600  ft.,  400  ft.  and 
300  ft.  from  the  tail  end  with  1850  ft.  of  cable  towed;  1400  ft.,  200  ft  and 
100  ft.  from  the  tail  end  with  1650  ft.  of  cable  towed;  and  400  ft.  and  300  ft. 
from  the  tail  end  with  1150  ft.  of  cable  towed.  It  can  be  seen  that  all  spec- 
tra lie  very  close  to  one  another.  For  practical  purposes,  therefore,  they 
can  be  represented  by  a single  spectrum. 

This  suggests  that  the  longitudinal  strain  energy  in  the  cable  is  supplied 
locally,  with  the  cable1 s own  turbulent  boundary  layer  being  the  prime  candi- 
date for  the  source,  rather  than  a vibrating  tow  point  or  the  strong  turbu- 


23 


LONGITUDINAL  STRAIN  (dBre  lm/mA/RF) 


Fig. 

100, 


15.  Longitudinal  strain  spectra  at  12  knots  at  gage-pairs  located 
200,  300,  400,  1400  and  1600  feet  from  the  tail  end. 
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lence  in  the  immediate  vicinity  of  the  submarine.  This  point  of  view  is  very 
well  corroborated  by  the  detailed  structure  of  the  spectrum  at  the  forward 
gage  pair.  Two  such  spectra  are  cluttered  up  with  six  others  in  Fig.  15,  but 
one  of  them  is  shown  alone  in  Fig.  16. 

The  structure  of  interest  here  is  the  sequence  of  regularly  spaced  maxima 
and  minima  in  the  spectrum.  By  projecting  the  sequence  back  to  the  origin, 
one  finds  that  a maximum  would  occur  there.  This  fits  a standing  wave  model 
in  which  the  strain  is  proportional  to  sin(knz),  where  z is  the  distance  from 
the  tow  point  and  k^  is  the  longitudinal  strain  wave  number.  Since  k^  = 0)/cf  , 
where  cr  is  the  longitudinal  strain  wave  velocity,  the  model  predicts  that  a 
strain  maximum  occurs  when  oo  = 0. 

The  model  also  predicts  that  the  minima  should  have  a constant  spacing 
Af  given  by  27TAfz/c^  = 7T.  Thus,  from  the  spacing  of  about  26.5  Hz  shown  in 
Fig.  16,  and  from  the  value  of  76  m (250  ft.)  for  z,  the  wave  speed  c^  works 
out  to  be  2170  m/sec.  But  c^  can  also  be  calculated  from  the  cable  longitud- 
inal stiffness  EA  and  its  effective  mass  per  unit  length  m^ . Then 
c,  = v^A/m^.  The  measured  value  of  EA  (see  Section  IV)  is  0.9  x 10^N  and  the 
cable  specific  gravity  is  about  0.75,  so  c^  turns  out  to  be  2050  m/ sec  by 
this  method.  The  close  agreement  between  the  two  speed  calculations  leads  one 
to  conclude  that  the  regular  structure  of  the  spectrum  arises  from  the  wave 
interference,  at  the  location  of  the  gage  pair,  between  the  strain  energy 
propagating  along  the  cable  towards  the  tow  point  and  the  same  energy  being 
reflected  at  the  tow  point  and  propagating  back  the  other  way. 

If  the  source  of  the  strain  energy  on  the  cable  were  a vibrating  tow 
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Fig.  16.  Longitudinal  strain  spectrum  at  12  knots  at  the  forward  gage-pair 
(1600  ft.  from  the  tail  end,  250  ft.  from  the  tow  point). 
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point,  the  energy  would  tend  to  flow  only  in  one  direction,  towards  the  tail. 
The  effect  of  reflections  at  the  tail  would  not  be  readily  observable  at  the 
forward  gage  pair  because  of  the  attenuation  involved  in  traveling  practically 
the  whole  length  of  the  cable  in  both  directions.  This  attenuation,  at  45  Hz 
for  example,  is  about  20  dB.  Thus  the  interference  effect  would  be  absent. 

The  characteristics  discussed  above  and  depicted  in  Figs.  15  and  16  apply 
equally  well  at  the  other  speeds  used.  That  is,  that  the  strain  spectrum  is 
independent  of  position  along  the  cable  and  derives  its  energy  of  excitation 
from  the  cable fs  own  turbulent  boundary  layer.  By  drawing  a single  smoothed 
average  spectrum  through  the  eight  spectra  shown  in  Fig.  15,  and  repeating 
this  for  the  other  speeds  the  smoothed  average  curves  shown  in  Fig.  17  were 
obtained.  These  show  the  spectral  density  of  the  longitudinal  cable  strain 
as  a function  of  speed  and  frequency. 

It  should  be  noted  that  these  strain  spectra  have  been  derived  from  the 
strain  gage  signals  using  a measured  value  for  the  gage  calibration.  The 
strain  spectrum  plotted  in  Ref.  2,  and  used  there  to  calculate  the  magneto- 
strictive  noise  voltage  spectrum  in  a certain  antenna,  was  based  on  a calcu- 
lated gage  calibration.  The  measured  calibration  was  obtained  by  mechanically 
exciting  the  section  of  cable  including  each  gage  pair  with  a known  alternat- 
ing longitudinal  strain  and  measuring  the  resulting  gage  amplifier  output 
voltages.  The  calculated  calibration,  on  the  other  hand,  simply  assumed  that 
the  gage  behaved  like  the  ideally  limp  gage  which  strains  the  same  amount  as 
the  surrounding  medium  without  constraining  it  in  any  way.  The  measured  cal- 
ibration resulted  in  a value  of  8 dB  re . lV/unit  strain,  whereas  the  calcu- 
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Fig.  17.  Smoothed  average  longitudinal  strain  spectra  at  various  speeds. 

The  longitudinal  strain  spectrum  is  essentially  independent  of  position  along 
the  cable  except  for  an  interference  effect  observed  near  each  end. 
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lated  calibration  was  14  dB  re.  lV/unit  strain.  That  is,  the  gage,  together 
with  its  mount,  is  sufficiently  stiff,  compared  with  the  soft  polymeric  jacket 
in  which  it  is  imbedded,  for  the  gage  strain  to  be  only  half  that  of  the  cable 
itself.  This  means  that  the  strain  spectrum  in  Fig.  13  of  Ref.  2,  and  the 
corresponding  voltage  spectrum,  should  be  raised  6 dB . 

If  the  longitudinal  strain,  with  its  long  wavelength  and  low  attenuation 
per  unit  length,  is  explained  by  a local  generation  model,  then  the  transverse 
vibration,  with  its  short  wavelength  and  high  attenuation  per  unit  length, 
would  be  expected  certainly  to  be.  The  results  of  the  test  support  this 
expectation.  Figure  18  shows  the  curvature  spectra  at  five  different  loca- 
tions along  the  cable.  (The  spectrum  for  the  300  ft.  location  is  omitted 
because  it  lies  more  or  less  directly  on  top  of  that  for  200  ft.)  At  high 
frequencies,  where  the  cable's  wave  properties  are  determined  by  the  cable 
stiffness  rather  than  its  tension,  the  spectra  all  merge  into  one.  At  low 
frequencies  this  does  not  occur,  because  there  the  tension,  which  changes 
continuously  over  the  length  of  the  cable,  becomes  predominant.  Whether  the 
local  generation  model  fits  well  at  the  lower  frequencies  can  be  examined  only  by 
using  the  cable's  mechanical  properties  (at  the  location  of  each  gage  pair 
and  at  each  speed)  to  deduce  the  driving  force  spectrum  S^(0,O))  from  the  cable 
response  (the  curvature  spectra) . This  will  be  done  in  the  next  section  for 
all  the  testing  speeds. 

Before  leaving  this  section,  however, it  is  necessary  to  discuss  the 
results  of  the  circular  turn  test  runs,  the  slow  speed-increase  test  run  and 
the  two-wire  gage  connection  test. 
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Fig.  18.  Cable  curvature  spectra  at  five  different  gage-pair  locations. 
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The  spectra  of  longitudinal  strain  and  of  curvature  measured  during  cir- 
cular turns  are  shown  in  Figs.  19  and  20.  If  they  are  to  be  believed,  they 
show  that  at  the  forward  gage  pair  (i.e.,  near  the  submarine)  both  kinds  of 
vibration  can  be  some  40  dB  or  more  higher  in  level  that  during  a straight  run 
They  are  believable  in  the  sense  that  this  is  true  for  at  least  one  of  the 
kinds  of  vibration.  But  since  the  data  reduction  procedure  depends  upon  long- 
itudinal strain  being  statistically  independent  of  curvature,  it  may  not  be 
true  for  both.  For  a curved  cable,  the  two  are  not  independent  and  so  the 
data  reduction  is  suspect.  It  would  not  grossly  magnify  the  apparent  energy 
of  both,  but  it  could  misapportion  the  energy  between  the  two. 

Near  the  tail  end  of  the  cable  (400  ft.  from  the  end,  to  be  precise),  the 
effect  of  the  circular  turn  is  not  nearly  as  marked  as  it  is  near  the  tow 
point.  In  fact,  Figs.  19  and  20  seem  to  show  that,  at  the  higher  frequencies, 
the  vibration  levels  are  sometimes  slightly  less  during  circular  turns  than 
during  straight  runs. 

The  result  of  the  circular  turn  tests  is,  therefore,  that  at  low  fre- 
quencies and  near  the  tow  point,  a circular  turn  induces  much  higher  vibration 
levels  In  the  cable  than  does  a straight  run  at  the  same  speed.  Whether  the 
vibration  spectra  during  turns  are  precisely  as  shown  in  Figs.  19  and  20  is 
in  doubt,  however,  because  the  data  reduction  program  may  have  misapportioned 
the  total  strain  energy  between  the  two  kinds  of  vibration.  . 

The  purpose  of  the  slow  speed-increase  test  was  to  apply  a slowly  increas 
ing  tension  to  the  cable  and  thereby  cause  the  helically  wound  strength  mem- 
bers on  the  cable  slowly  to  unwind.  This  would  then  cause  the  cable  at  each 
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Fig.  19.  Longitudinal  strain  spectra  at  different  speeds  during  circular 
turns . 
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Fig.  20.  Transverse  curvature  spectra  at  different  speeds  during  circular 
turns . 
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gage  location  to  revolve  slowly  about  its  longitudinal  axis.  Then,  if  at  one 
moment  a gage  pair  were  registering  the  cable  curvature  in  a vertical  plane, 
some  time  later  they  would  be  registering  the  curvature  in  the  horizontal 
plane.  In  this  way,  an  azimuthal  asymmetry  in  the  transverse  vibration  of  the 
cable  would  be  observable  as  a cyclic  variation  in  the  mean  square  output 
voltage  of  the  gages. 

In  fact,  no  such  cyclic  variation  was  observed.  This  could  be  explained 
by  the  proposition  that  the  gage  signal  arising  from  cable  curvature  was  much 
less  than  that  arising  from  longitudinal  cable  strain,  and  so  the  curvature 
signal  was  unobservable  to  begin  with.  In  fact,  the  reverse  was  the  case. 

The  gage  signal  arising  from  the  cable  curvature  was  at  least  13  dB  greater 
than  that  arising  from  the  longitudinal  strain.  Thus  it  appears  that  the 
transverse  vibration  of  the  cable  is  statistically  independent  of  the  plane  in 
which  it  is  measured. 

The  last  test  to  be  carried  out  was  an  experiment  to  see  if  the  three-wire 
gage  connection  used  for  all  data  runs  was  really  necessary.  It  was  used  as 
a precaution  against  the  possibility  of  the  strain  in  the  long  connecting 
wires  producing  a spurious  strain  signal  [6].  Figures  21  and  22  show  the 
vibration  spectra  obtained  using  a two-wire  gage  connection  with  those 
obtained  using  the  normal  three-wire  connection.  The  necessity  for  using  the 
latter  is  clear.  The  error  incurred  in  using  a two-wire  connection  is  seen 
to  be  more  for  the  common  mode  longitudinal  strain  signal  than  for  the  dif- 
ferential mode  curvature  signal.  This  is  to  be  expected,  for  the  cable  con- 
ductors all  strain  in  unison  and  therefore  tend  to  be  self  cancelling  in  the 
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Fig.  21.  Longitudinal  strain  spectrum  measurement  using  a two-wire  con- 
nection compared  with  that  using  the  normal  three-wire  connection. 
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Fig.  22.  Transverse  curvature  spectrum  measurement  using  a two-wire  gage 
connection  compared  with  that  using  the  normal  three-wire  connection. 
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error  they  cause  in  the  strain  difference. 

The  comments  made  earlier  concerning  the  measured  versus  calculated  gage 
calibration  for  longitudinal  strain  apply  also  to  the  calibration  of  the 
gages  to  cable  curvature.  The  measured  curvature  calibration  is  -40  dB 
re.  lV/m  \ which  is  6 dB  less  than  the  calculated  calibration.  Thus  the 
spectra  of  curvature  and  derived  mo tion-induced  noise  appearing  in  Fig.  13  of 
Ref.  2 should  be  raised  by  6 dB . 

IV . Derived  Force  Spectra 

The  smoothed  average  longitudinal  strain  spectra  of  Fig.  17  are  seen  to 

-2  -4 

vary  approximately  as  f and  v , where  v is  the  towing  speed.  But  if 
Sg(0,u)  is  the  two  dimensional  power  spectral  density  of  the  longitudinal 
shear  force  per  unit  length  on  the  cable,  then  according  to  the  local  exci- 
tation model  of  cable  vibration  [2],  the  strain  spectrum  (to)  is  given  by 

, \\ 

Sx(u)  = k^Sx(w)  = -------  Ss(0,co)  (1) 

TTD  |m  r.  I 

1 X,  1 

Here  is  the  Q-factor  of  the  cable  vibration,  k - is  the  longitudinal  strain 

wave  number  and  ni,  the  effective  mass  per  unit  length  for  longitudinal  vibra- 

x, 

tion.  All  these  parameters  are  defined  in  Ref.  2. 

Since  Q.  is  a weak  function  of  frequency,  and  since  is  proportional 

_2 

to  frequency,  then  to  arrive  at  the  result  that  (lj)  is  proportional  to  f “ , 

it  must  be  the  case  that  S^CO,^)  is  proportional  to  f Therefore,  to  find 

a quantitative  expression  for  the  longitudinal  shear  force  spectrum  S (0,uj)  which 

s 

fits  the  measured  cable  strain  it  is  appropriate  to  set 
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(2) 


4 

Sg (0,aj)  = a ~ 

and  choose  a so  that  the  resulting  strain,  calculated  using  (1),  is  a good 

fit  to  the  measured  curves  (Fig.  17). 

To  carry  out  the  parameter  fit  it  is  necessary  to  know  EA,  the  cable’s 

longitudinal  stiffness,  in  magnitude  and  phase,  and  also  the  specific  density 

b of  the  cable.  These  were  measured  to  be 

EA  = 0.9  x 106  (1  - iO . 025)  Newtons 
b = 0.75. 

-Q 

With  a value  of  a of  9.0  x 10  MKS  units,  the  strain  spectra  derived 
from  the  longitudinal  shear  force  spectrum  of  Eq.  (2)  are  shown  plotted  as 
solid  lines  in  Fig.  23.  The  dashed  lines  are  the  experimentally  measured 
smoothed  average  spectra  shown  in  Fig.  17.  It  can  be  seen  that  the  fit  is 

reasonably  good  for  practical  purposes. 

Thus,  given  the  forcing  function  spectrum  (2),  one  can  calculate  the 

resulting  cable  strain  S^(a))  for  any  0.65  inch  diameter  cable  by  simply  mea- 
suring EA  and  b for  the  cable  in  question.  The  only  other  condition  necessary 
for  the  validity  of  this  derivation  is  that  the  surface  of  the  cable  should 
be  hydrodynamically  smooth  [2] — a fairly  benign  requirement  for  the  Reynolds 
numbers  encountered  in  this  application. 

2 

There  is  good  reason  to  believe  [7]  that  Sg(0,a3)  is  proportional  to  a . 

If  that  is  so,  then  the  formula  for  it  can  be  generalized  further  to 

,42 

S (0,u>)  = 1.32  X 10  , (3) 

all  units  MKS.  Then  the  strain  spectrum  in  a cable  of  arbitrary  diameter  and 
mechanical  properties  can  be  calculated. 
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Fig.  23.  Strain  spectra  calculated  using  the  simple  longitudinal  shear 
force  spectrum  given  in  Eqs.  (2)  or  (3)  (solid  curves)  compared  with  smoothed 
average  experimental  curves  (dashed  curves) . 
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For  transverse  cable  vibration,  the  local  excitation  model  [2]  predicts 
that  the  curvature  spectrum  S ^ (co)  is  given  by 

A k 5Q 

S (u>)  = kt  S (u>)  = j sf(°»u)-  (4) 

y TTU)  | m I 

Here  k^_,  and  m^_  are  the  wave  number,  Q-factor  and  effective  vibrating  mass 
associated  with  transverse  cable  vibration  and  S^(0,oj)  is  the  two  dimensional 
power  spectral  density  of  the  transverse  force  per  unit  length  on  the  cable. 
From  this  equation,  if  Sc(oj)  is  known,  then  S^(0,oo)  can  be  calculated.  The 
only  extra  data  necessary  are  the  buoyancy  b and  cable  bending  stiffness  El. 

From  these,  using  the  formulas  of  [2],  k^_,  Q and  m^  can  be  calculated. 

The  cable  bending  stiffness  El  was  measured  to  be 

El  = 1.2(1  iO.12)  Newton  m^ 

and  the  buoyancy,  as  mentioned  before,  was  measured  to  be  0.75.  Using  these 
values,  the  transverse  force  spectrum  S^(0,o))  corresponding  to  each  measured 
curvature  spectrum  (oo)  was  calcuated.  The  result  is  the  lightly  drawn 
sets  of  spectra  in  Fig.  24. 

There  is  a disturbingly  large  scatter  in  Fig.  24  in  the  different  spec- 
tral density  levels  calculated  for  a particular  frequency  and  speed.  However, 
the  derivation  of  S^(0,oj)  from  Sc(oo)  involved  what  may  be  an  over  simplified 
vibration  model  in  which  the  bending  stiffness  is  independent  of  tension  and 
frequency  and  which  requires  that  the  tension  be  an  accurately  known  function 
of  speed. 

Another,  more  direct,  method  of  deriving  S^(0,oo)  is  by  way  of  the  elec- 
trode-pair antenna’s  measured  noise  voltage.  If  the  noise  voltage  is  motion- 
induced,  then  it  is  related  to  S^(0,oo)  by  the  following  simple  law  [2] 
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Fig.  24.  Transverse  force  spectrum  derived  from  strain  gage  curvature  mea- 
surements (light  curves)  and  from  electrode  pair  noise  voltage  measurements 
(heavy  curves) . 
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(5) 


s (u)  = n 


sf(0,u). 


Here  2H  is  the  length  of  the  antenna,  B^_  the  transverse  component  of  the  geo- 
magnetic field  and  m^  the  effective  mass  of  the  cable  per  unit  length  for  trans- 
verse vibration.  But  the  measured  noise  voltage  in  a 0.65  inch  diameter, 

A 

300  m long  electrode-pair  antenna  in  a region  where  % 0,5  x 10  T is  very 
well  represented  as  a function  of  frequency  and  speed  by  [3] 


S 

v 


(03)  = 9.69  x 10 


-14 


(6) 


all  units  MKS . Thus,  assuming  the  antenna  cable’s  buoyancy  was  0.75,  one  can 
invert  (5)  and  (6)  to  obtain 


Sf (0,0))  = 7.16  x 10 


(7) 


all  units  MKS. 

This  spectrum  is  also  plotted  in  Fig.  24  as  the  five  heavily  drawn  curves. 

Comparing  the  sets  of  lightly  drawn  curves  in  Fig.  24  with  those  heavily 
drawn,  one  sees  that  the  total  range  as  a function  of  speed  is  about  the  same 
for  both  sets  and  the  general  trend  with  frequency  is  much  the  same.  Since 
the  agreement  is  so  marked,  it  must  be  the  case  that  the  electrode-pair  ant- 
enna’s dominant  noise  voltage  is  motion-induced  and,  since  the  law  connecting 
its  noise  voltage  with  the  transverse  force  spectrum  is  so  simple,  that  the 
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force  spectrum  can  be  accurately  derived  from  the  electrode-pair  antenna1  s 
noise  voltage.  That  force  spectrum  is  given  by  Eq . (7).  It  gives  the  fitted 
curves  shown  in  Fig.  1,  Ref.  3 when  substituted  in  (5). 

The  deviations  in  Fig.  24  of  the  force  spectra  derived  from  the  curvature 
measurements  from  the  force  spectra  derived  from  electrode-pair  measurements 
are  therefore  attributed  to  the  uncertainties  concerning  cable  tension,  stif- 
fness and  damping  as  functions  of  speed  and  frequency.  These  are  needed  to 

derive  force  from  curvature,  but  not  to  derive  force  from  voltage. 

2 

Again,  since  S^(0,(jo)  is,  like  S ^ ( 0 , to) , proportional  [7]  to  a , the  force 
spectrum  given  by  (7)  can  be  generalized  to 

2 5 2 

S (0,oj)  = 1.05  X 10  — - . (8) 

f 

From  this, and  the  properties  of  the  cable,  the  curvature  spectrum  of  an 
arbitrary  cable  towed  at  an  arbitrary  speed  can  be  calculated. 

V.  Conclusions 

The  accelerometer  and  strain  gage  measurements  have  shown  that  the  vibra- 
tion of  a cable  towed  behind  a submarine  is  locally  excited.  That  is,  the 
vibration  is  excited  by  the  fluctuating  shear  and  pressure  forces  on  the  sur- 
face of  the  cable  arising  from  the  turbulent  water  flow  in  the  boundary  layer 
surrounding  the  cable.  By  fitting  simple  expressions  for  the  shear  and  pres- 
sure force  spectra  to  the  measured  data,  quantitative  formulas  have  been 
derived  from  which  one  can  calculate  the  spectra  of  longitudinal  strain  and 
curvature  for  an  arbitrary  cable  towed  at  an  arbitrary  speed.  The  only  extra 
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data  one  needs  are  the  mechanical  properties  and  dimensions  of  the  cable. 

The  answers  to  the  questions  raised  in  the  Introduction  are  therefore  that 

i)  The  tail  need  be  no  longer  than  100  feet,  possibly 
could  be  even  smaller.  No  vibration  measurements  were  made  closer 
than  100  ft.  from  the  end  of  the  cable,  and  so  how  closely  within 
the  last  100  ft.  the  statistically  uniform  vibration  measured  at 
other  points  approaches  the  tail  end  is  not  yet  known. 

ii)  Since  both  longitudinal  strain  and  curvature  are  locally 
generated,  isolators  would  have  marginal  effect  on  the  cable  vibra- 
tion in  the  antenna  region. 

iii)  Since  the  cable fs  own  turbulent  boundary  layer  is  the  dom- 
inant source  of  the  cable  vibration,  moving  the  antenna  region  fur- 
ther away  from  the  submarine  wake  would  have  no  effect. 

iv)  The  quantative  agreement  between  the  measured  noise  voltage 
of  an  electrode-pair  antenna  and  its  noise  voltage  calculated  from 
the  vibration  measurements,  strongly  suggests  that  vibration,  rather 
than  hydromagnetic  effects,  is  the  source  of  the  noise. 
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APPENDIX  A 


STRAIN-GAGE  PAIR  DATA  PROCESSING 

If  the  amplifier  gains  were  precisely  their  nominal  values,  and  if  the 
same  were  true  for  the  gage  sensitivities  to  curvature  and  to  longitudinal 
strain,  then  the  output  voltages  v^  and  of  the  two  balanced  gage  channels 
would  be  given  by 


V1  = VA  + vc 
v2  = VA  “ Vc 


(Al) 


where  and  are  the  output  voltages  due  to  longitudinal  strain  and  curva- 
ture separately.  The  gains  may  not  be  their  nominal  values,  however,  and  the 
gage  sensitivities  can  also  be  different.  Thus  the  true  situation  can  be 
represented  by  replacing  Equations  (Al)  by 

101 

v,  = (av,  + 8v  )e 
1 Ac 

102 

v2  = (6vx  - Yvc)e  , 

where  a,  8,  6 and  y are  all  real  numbers  close  to  unity  to  account  both  for 
the  non-nominal  gage  sensitivities  and  for  the  non-nominal  gain  amplitudes, 
and  6^  and  are  real  numbers  to  account  for  the  non-nominal  gain  phases. 

Since  the  sensitivity  of  each  gage  to  cable  curvature  can  vary  independ- 
ently of  its  sensitivity  to  longitudinal  cable  strain,  8/a  is  not  necessarily 
equal  to  y/6.  Thus  all  four  independent  amplitude  factors  a,  8,  Y and  6 are 

needed.  The  phase  response  of  the  amplifier  affects  both  v^>  and  v equally, 

A c 
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however,  and  so  only  two  independent  phase  factors  0^  and  0^  are  needed. 

If  a,  6,  6,  y,  0^  and  0^  were  known  precisely,  and  could  be  eval- 
uated separately  from  any  given  pair  of  values  of  v^  and  v^  by  inverting  (A2) . 
Careful  measurement  could  fairly  readily  establish  0^  and  0^  with  sufficient 
accuracy,  but  an  accurate  measurement  of  the  amplitude  factors  would  require 
accurately  known  combinations  of  longitudinal  strain  and  curvature  to  be  estab- 
lished in  the  cable  while  each  gage  amplifier^  output  voltage  was  simultan- 
eously being  recorded.  The  mechanical  problems  involved  in  such  a venture 
are  formidable.  There  is  incentive  to  seek  an  easier  way. 

Thus  the  measurement  problem  consists  of  estimating  the  separate  power 

spectral  densities  and  S of  v.  and  v when  only  v.  and  v~  are  directly 

Aaccac  y 1 2 

available,  when  the  four  amplitude  factors  are  known  only  to  the  extent  that 
they  are  close  to  unity  and  when  the  two  phase  factors  are  known  to  be  small 

compared  to  a radian.  The  one  crucial  additional  fact  is  that  v.  and  v are 

A C 

statistically  independent  of  one  another — S.  = 0.  This  allows  the  separa- 

Ac 

tion  of  S,.  and  S to  be  carried  out,  even  though  one  might  be  much  larger 
A A CC 

than  the  other. 

Since  the  amplitude  factors  are  close  to  unity,  a measurement  of  the 

geometric  mean  quantities  /a6*v^and  /67vc  is  essentially  equivalent  to 

a measurement  of  v.  and  v . Thus  if  these  modified  voltages  be  represented 

as  v\  and  v?  , and  then  substituted  into  (A2)  , the  result  is 
A c 

i0i 

v = (av'i  + bv'  )e 
1 Ac 


v 


2 


1 

b 


v 


(A3) 
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where  a = /a/6‘  and  b = /g/y'. 

If  b were  equal  to  a,  then  since  v'.  and  v'  are  statistically  independ- 

A C 

ent,  the  spectral  density  of  would  be  times  the  spectral  density  of 
^22  V2  * ^ Edition,  an<^  were  zero>  the  separate  measurement  of 

and  would  then  involve  first  computing  the  spectral  densities  and 


cc 


S22  °f  V1  and  v^,  evaluating  the  amplitude  correction  h by  means  of  the  form- 
ula h^  = aUC^  t^en  computing  the  spectral  densities  of 

[(l/h)v1  + hv^l/2.  These  last  are  equal  to  and  S , respectively. 

A procedure  very  similar  to  this  can  be  used  when  b £ a and  when  0^, 

&2  ^ 0.  The  amplitude  correction  is  calculated  in  the  same  way  (and  is 
therefore  a compromise  between  the  unequal  correction  factors  needed  to  deal 
with  the  unequal  amplitude  factors  a and  b)  and  the  phase  correction  is 

obtained  by  using  the  cross  spectral  density  (which  is  the  expected  value 

* 

of  v^v2  The  derivation  of  the  procedure  is  as  follows: 

The  correction  factor  h is  given  by 
S, 


4 _ “11 
h - s 

22 


(A4) 


and,  ideally,  the  phase  differences  between  the  channels  should  be  corrected 

-i0  -i02 

by  modifying  v^  and  v^  by  e and  e . That  is, one  computes  the  spectral 


densities 


-ie  -ie9 

s+  ■ <l^h  Vle  - hv2e  ^ > 


(A5) 


where  the  angle  brackets  denote  expected  value.  (It  will  transpire  that  0.. 
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and  need  not,  in  practice,  be  evaluated  directly.  See  below.) 

When  (A4)  is  substituted  in  (A5)  and  the  expected  value  of  the  sum  of 
terms  rewritten  as  the  sum  of  the  expected  values  of  the  terms,  the  result 
is 


S+  'I  l/S7S  ± <SU-Scc>'- 


(A6) 


But  since  = ^vlv2  / = )e  > an d since  0^  and  are 

assumed  to  be  small  compared  with  a radian,  the  term  - S is  equal  to 
sgn(Re{S^2^)  Is^U  Therefore,  (A6)  can  be  rewritten  as 


S+  = 2 [/SllS22'  ± 8gn(Re{S12))|S12|]. 


(A7) 


All  the  quantities  in  this  expression  can  be  computed  from  v^  and  v^ . 
and  S^2  are  the  auto  spectra  of  v^  and  and  S.^  is  the  cross  spectrum. 

The  question  remaining  is  the  extent  to  which  S+  and  S as  defined  by  (A7) 

are  estimates  of  and  S 

AA  cc 

By  direct  algebraic  manipulation  of  (A3),  it  can  be  shown  that 


/S11S22  ' (SXX  + Scc>  {1  + 


s,,s 

AA  CC 


1/2 


(SXX  + Scc>' 


(A8) 


where  £ is  the  "asymmetry  factor" 


5 = 


l 

4 


(A9) 
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Therefore,  by  combining  (A6)  and  (A8) , and  expressing  the  square  root  as  the 
first  two  terms  of  its  Taylor  series,  one  obtains 


cc 


S+  SAA(1  + S, ,+S 

AA  cc 


5 + o{C» 


s = s (i  + 

_ CC  S-.-.+S 

AA  CC 


xx  s + oU2}). 


(A10) 


These  expressions  show  that  the  sum  and  difference  spectra  and  S , evalu- 
ated by  means  of  the  formula  (A7) , are  closely  equal  to  the  required  strain 

and  curvature  spectra  and  S . This  remains  true  for  both  of  them  even 

AA  cc 

when  one  is  very  large  compared  with  the  other.  In  fact,  by  examining  the 
exact  formulas  for  which  (A10)  are  an  approximation,  one  can  show  that  in 
general 


<1+C>su 


s 

cc 


< S < (l+£) s 
— _ — cc 


(All) 


If  the  unlikely  event  should  occur  that  a/b  is  as  far  from  its  nominal  value 
of  unity  as  1.5  (or  0.667),  then  C = 0.17,  and  in  that  case  the  maximum  pos- 
sible error  in  using  S,  and  S as  estimates  of  S, , and  S is  17  percent  or 
0.68  dB. 

An  example  of  the  process  is  illustrated  in  Figs.  A1  through  A6.  The 
first  four  show  the  two  auto  spectra  S^(f)  and  $22(0  followed  by  the  cross 
spectrum  S^CO  in  amplitude  and  phase.  These  spectra  were  taken  using  1024 
samples  per  second  of  the  analog  strain  signals  from  gages  3 and  4 during  a 
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16  knot  straight  run.  The  number  of  samples  used  for  each  Fourier  transform 
was  512,  which  implies  a resolution  of  2 Hz.  The  mean  square  of  a total  of 
128  such  transforms  comprized  each  auto  spectrum.  The  cross  spectrum  was  the 
mean  of  128  cross  products.  As  large  a number  as  128  transforms  were  aver- 
aged in  order  to  get  good  estimates  of  S-^(f)  , S22^  an d Only  then 

could  a small  longitudinal  strain  spectrum  S^(f)  be  accurately  recovered  in 
the  presence  of  a large  curvature  signal  Scc(f). 

The  result  is  shown  in  Figs.  A5  and  A6,  which  give,  respectively,  the 
sum  spectrum  S_^_(f)  and  difference  spectrum  S (f). 

Notable  features  present  in  the  figures  are: 

i)  The  phase  between  the  two  strain  signals  (Fig.  A4) . It 
stands  very  close  to  180°,  except  where  perturbed  by  line  fre- 
quency interference,  across  the  whole  band  of  frequencies.  This 
indicates  that  the  strain  signals  are  dominated  by  the  curvature 
signal,  which,  being  a difference  mode  signal,  is  180°  removed 
in  phase  in  one  gage  with  respect  to  the  other. 

ii)  The  sum  spectrum  S+(f)  (Fig.  A5) , which  is  the  estimate 
of  the  longitudinal  strain  spectrum  S^^(f),  has  a much  lower  level 
than  the  estimate  S (f)  of  the  curvature  spectrum  Scc(f)  (Fig.  A6). 

This  is  consistent  with  the  phase  data  of  Fig.  A4 . Also,  the  sum 
spectrum  S+(f)  shows  strong  line  frequency  interference  whereas  the 
difference  spectrum  S (f)  does  not.  This  would  be  expected,  since 
line  frequency  interference  is  a common-mode  phenomenon  for  iden- 
tical gage  connections  to  identical  amplifiers.  Particularly 
gratifying  about  the  two  figures,  however,  is  that  the  separation 
between  the  two  signal  sources  is  so  consistent.  In  places,  the 
spectra  differ  by  more  than  two  orders  of  magnitude,  and  the  comp- 
utation of  spectral  levels  is  computed  anew  at  every  increment  of 
2 Hz.  Nevertheless,  the  resulting  spectra  are  smooth  and  distinct. 

iii)  The  cross-spectrum  amplitude.  Fig.  A3,  is  apparently 
indistinguishable  from  either  of  the  auto  spectra.  Figs.  A1  and 
A2,  and  yet  the  difference  between  them,  in  the  form  expressed 
by  (A7) , is  not  simply  a random  error.  It  is  the  sum  spectrum 
S_(f)  plotted  in  Fig.  A5 . (Since  the  phase  of  Si2(f)  is  180°, 
its  real  part  has  a negative  sign.) 
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Fig.  Al.  The  auto-spectrum  of 
16  knots  and  steady  heading. 


the  strain  signal 


from  gage  3 at  a speed  of 
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Fig.  A2 . The  auto-spectrum  of  the  strain  signal  from  gage  4 at  a speed  of 
16  knots  and  steady  heading. 
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Fig.  A3.  The  amplitude  of  the  cross-spectrum  of  the  strain  signals  from 
gages  3 and  4 at  a speed  of  16  knots  and  steady  heading. 
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Fig.  A4.  The  phase  of  the  cross-spectrum  of  the  strain  signals  from  gages 
3 and  4 at  a speed  of  16  knots  and  steady  heading. 
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Fig.  A5.  The  estimate  of  the  longitudinal  strain  spectrum  at  the 
of  gages  3 and  4 at  a speed  of  16  knots  and  steady  heading. 
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Fig.  A6.  The  estimate  of  the  curvature  spectrum  at  the  location  of  gages 
3 and  4 at  a speed  of  16  knots  and  steady  heading. 
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